Ergodicity breaking in a model showing many-body localization 
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We study the breaking of ergodicity measured in terms of return probability in the evolution of a quantum 
state of a spin chain. In the non ergodic phase a quantum state evolves in a much smaller fraction of the Hilbert 
space than would be allowed by the conservation of extensive observables. By the anomalous scaling of the 
participation ratios with system size we are led to consider the distribution of the wave function coefficients, a 
standard observable in modern studies of Anderson localization. We finally present a criterion for the identifi- 
cation of the ergodicity breaking (many-body localization) transition based on these distributions which is quite 
robust and well suited for numerical investigations of a broad class of problems. 



The question of whether an Anderson localization (AL) 
transition [ I ] can occur in a system of interacting particles has 
been recently suggested to have a positive answer [2, 3]. The 
mechanism which underpins this effect (dubbed many-body 
localization or MBL transition) requires the interaction to act 
in a substantially non-perturbative way, therefore providing an 
example of how disorder and strong interactions interplay in a 
quantum theory. 

The natural setup to study the MBL transition is the dy- 
namics (these were also the terms of the question posed in 
[" ]) and in this perspective it is a question about the founda- 
tions of statistical mechanics, namely, on the validity of the 
ergodic hypothesis. MBL also presents the terms in which a 
quantum glass can be defined and from there it is only a small 
leap to conjecturing that MBL is a natural ingredient for hard 
computational quantum problems [ , ] (as Ising spin glasses 
are a natural scenario to discuss the Physics of hard combina- 
torial optimization problems [6, 7]). Even by neglecting the 
implications for experiments (and there are many [8, 9]) the 
topic should be considered worth of serious investigation. 

One dimensional systems are particularly suited for study- 
ing MBL because the single particle spectrum is completely 
localized for arbitrarily small disorder, and therefore any ob- 
servation of derealization must be attributed to the interac- 
tion. In this paper, we analyze in detail the ergodicity prop- 
erties of an XXZ chain with random fields. This particu- 
lar example has already provided different indications of the 
MBL transition for sufficiently large disorder: in [10] corre- 
lation functions and spectral properties were studied, while in 
[11, 12] tDMRG was used to investigate the different satura- 
tion properties of the entanglement entropy in the two phases. 
While the existence of a transition in the dynamics of this 
model is now almost certain, its precise location, the possible 
existence of a critical phase and the nature of the phases that it 
separates are subject of debate. This should not be regarded as 
a debate about a particular spin chain but rather as an attempt 
at characterizing as much as possible the differences between 
MBL and AL. 

Consider the real time evolution of a state \i/j ) as it is en- 



coded into the Green's function (h = 1 in the rest of the paper) 

G(t) = (^ \e- ttH |Vo) ■ (1) 
We introduce the inverse participation ratios as the moments 

IPR, = ^| Wo)| 29 . (2) 

E 

where the sum runs over the full set of eigenstates \E). The 
long time average of the survival probability (the average re- 
moves some finite size effects like quasi-periodicity etc.) can 
be expressed as 

1 f T 

P= lim - / dt \G(t)\ 2 = IPR 2 . (3) 

Here (IPR 2 )~ 1 is therefore a measure of the portion of ex- 
plored Hilbert space during the quantum dynamics and it is 
usually dubbed participation ratio (PR). Analogously, higher 
order IPR 9 's describe finer details of the dynamics. 

Let us now comment on the choice of a suitable initial state 
for a Gedankenexperiment aimed at testing the breaking of 
ergodicity. First of all, consider what happens if we take 
a random state in the Hilbert space (therefore not an eigen- 
state) conditioned just to have an expectation value of the en- 
ergy E (with high probability, for a random state and a lo- 
cal Hamiltonian the standard deviation <5 = O (/Y 1 / 2 ) <C 
E = O {N)). The average values of local operators in this 
state will not show signs of ergodicity breaking. In fact, even 
at very large disorder there are states very close in energy 
(AE = O (e~ s ), where S is the microcanonical entropy at 
energy E) which are macroscopically different and the expec- 
tation value of a local operator will be the average of its values 
in these localized eigenstates, concealing the effect of disor- 
der (as expected from the ergodic theorem [13]). If we want 
to observe the effect of disorder on the dynamics, a reasonable 
prescription consists in choosing an eigenstate of the part of 
the Hamiltonian which dominates in the strong disorder limit. 
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Starting the dynamics coincides then with turning on the rest 
of the Hamiltonian. In the delocalized phase, during the quan- 
tum dynamics, the motion covers a finite fraction of the full 
Hilbert space (each eigenstate being individually thermal, the 
so-called "eigenstate thermalization hypothesis" (ETH) [14- 
1 6]). Instead, in presence of strong disorder, ergodicity breaks 
down and the many-body wave function motion is constrained 
on a small section of the full Hilbert space. 

We also believe that this point of view on MBL is what bet- 
ter brings forward its implications for quantum computation 
(or at least for the performance of the Adiabatic Algorithm 
[17]). In the localized phase the system gets frozen, the dy- 
namics unable to efficiently explore the Hilbert space, so the 
algorithm is not efficient in finding the ground state [4, 5, 17]. 

This view on the MBL transition will be the focus of this pa- 
per. We will show how the usual criteria for detecting AL need 
to be tweaked to capture the MBL transition; we will study the 
IPR's and will show how, although much information is con- 
tained in them, it is actually necessary to study the distribution 
of wave-function coefficients (ijj^\E), which is heavily tailed 
both in the localized and delocalized regions. 

We consider the Hamiltonian 
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with periodic boundary conditions. As the Hamiltonian com- 
mutes with the total z spin S z = ^ sf, we focus on the 
subspace with S z = 0. The random fields are chosen from a 
box distribution hi € [—h, h\. The model can be cast into a 
theory of fermions {S z = corresponds to half-filling), with 
on-site disorder hi. 

The As z s z term can be written as a two-body, point-like 
interaction for the fermions and for zero temperature it can be 
included perturbatively or non-perturbatively [ ] leading to 
an interesting phase diagram. When A = the fermions are 
free, an arbitrarily small disorder localizes the entire spectrum 
and therefore ergodicity is broken for any h > h c = 0. As 
A is increased MBL would appear as a peculiar phase tran- 
sition (possibly even at infinite temperature) at a critical h c 
increasing away from zero. On the other hand, for A> J the 
disorder necessary to break ergodicity should decrease again. 
In fact, for large A the relevant degrees of freedom are the 
domain walls of the classical Ising chain obtained by setting 
J = in (4). Longer domain walls have smaller hopping 
matrix elements and therefore they are more prone to local- 
ization than the fermions at J A. Once a few of these 
large domain walls have frozen, ergodicity can be considered 
broken and this occurs for smaller h, since both the effective 
hopping and interaction are smaller (effective randomness is 
always h). Here we present results of exact diagonalization 
for A = J = 1, where the delocalized phase is largest. 



According to the discussion of the previous section, we 
should test ergodicity by taking an initial state tp as one of 
the jV = (^2) configuration of spins |a) polarized along the 
z or — z direction, (e.g. \a) = ...)). We need to stress 
a major difference in the behavior of IPR2 in the localized 
and delocalized phases between AL and the present situation. 
While in the former one can distinguish the two phases by the 
participation ratio being 0(1) or not in the thermodynamic 
limit, this is not a sufficient criterion for us. For a many -body 
state, even in absence of interaction, IPR2 will be exponen- 
tially small in N also in presence of strong disorder, simply 
because each degree of freedom will have a localization length 
small but finite, corresponding to an individual participation 
ratio smaller than 1: multiplication of O (N) of these factors 
leads to an exponentially small IPR2. We need to correct the 
previous criterion by requiring that the delocalized and local- 
ized phase are distinguished by whether the ratio IPR2 /W _1 
is O (1) or not. The other IPR g 's, properly rescaled with pow- 
ers of the Hilbert space dimension Af, also represent indica- 
tors of ergodicity breaking. 

However, as far as averages over the initial states are in- 
volved we have found that PR's have better finite-size behav- 
iors (more on this later), so we considered: 



(5) 



{h},a 



where the subscripts in the average correspond to disorder 
realizations (indicated with h) and initial spin configuration 
|a)[19]. In particular the data for I2, shown in Fig. 1, are con- 
sistent with the lirriAr_j. 00 I^ih) = 12(h) where 12(h) = 0, 
for h > h c = 2.7 ± 0.3, although the finite-size corrections 
are strong already at h > 1.5. The prediction of 1/3 for 12 
at small h coming from the GOE ensemble, although qualita- 
tively correct, is quantitative inaccurate. This could be how- 
ever a consequence of the many-body structure in the finite 
size scaling that we did not take into account up to now. 
A similar information is obtained by the diagonal entropy 
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which is plotted for varying h in Figure 2 and also this quan- 
tity is clearly far from its thermodynamic limit of S = 1 in 
the delocalized phase. If we identify the critical point (see the 
arrows in Fig. 1 and Fig. 2) as the place where the N depen- 
dence sets in (for I 2 ) or drops out (for S) then both quantities 
identify a critical point consistent with h c = 2.7 ± 0.3 consis- 
tently with the findings of [10]. 

The diagonal entropy and the IPR's show that the wave 
function covers a number of sites that grows exponentially 
with the system size, although the exponent is smaller than in 
the ergodic phase. This suggests that in a many-body system, 
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We will consider therefore the average over eigenstates, initial 
states and disorder realizations: 



(x,N) = (6(x-Af\ (a\E)f)) a E ^ 



{h} 



(7) 



In the following we will drop the subscripts in the averages. 
This function depends both on x and N in general but in the 
ergodic delocalized phase, as Af plays the role of the space 
volume, we see that the dependence on N drops out [20-22] . 
We can then write the various IPR's as 



/>oo 

(IPR 9 ) =Af 1 - q / dx 
Jo 



x q (j){x) 
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FIG. 1. Average fraction of occupied Hilbert space as a function of h 
for different system sizes N — 8 to 16 using exact diagonalization. 
Notice how the limit for h — > is different from 1/3 which is the 
RMT prediction. 




FIG. 2. Average diagonal entropy as a function of disorder strength 
for different sizes N = 8, 10, 12, 14, 16. From the iV-dependence 
the transition is identified at h c ~ 2.7 ± 0.3. 
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the localized phase is necessarily characterized by the break- 
ing of ergodicity, but not necessarily by a concrete localiza- 
tion (IPR2 — O(l)). However, to pinpoint the transition and 
understand the reasons of this anomalous scalings we should 
analyze the full probability distribution of | (a\E) | 2 . 



DISTRIBUTION OF WAVE FUNCTION AMPLITUDES 

If one considers the various IPR g averaged over \a), one 
observes a peculiar scaling with N of each of them, which 
can be considered as due to large fractal dimensions. In this 
scenario, the safest observable to consider is the distribution 
of the properly rescaled wave function coefficients. As we are 
interested in typical states (infinite temperature) we will not 
follow the usual route of fixing the energy of the state but we 
will rather integrate over the whole spectrum. In the thermo- 
dynamic limit this corresponds to energy density E/N = 0. 



FIG. 3. The distribution <j> of scaled wave function amplitudes x = 
A/ - 1 (a\E) | 2 for different values of h. Upper panel: h = 1.2 in the 
middle of the ergodic phase where the scaling is perfectly verified, 
lower panel h — 4.2 in the many-body localized phase. In each 
figure the different curves correspond to different values of N, from 
8 to 16. Each curve is obtained by binning of not less than 3 10 6 
squared amplitudes. 



Illustrative plots are shown for different regimes in Fig. 3. 
As we said, even though in the ergodic phase, with this scaling 
the curves for different sizes collapse (similarly to AL), the 
distribution has an elbow at x ~ 1 and we find 



(x) oc 
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(9) 



where a, (3 depend on h. We have a < 1 < (3 ensuring the 
normalization of the distribution function in the delocalized 
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FIG. 4. The value of the exponent a (blue squares) and /3 (pink 
circles) in Eq. (9) for N = 16 (these exponents are independent 
of N within the symbol size). The exponent a crosses the value 1 
required by summability, which occurs at h ~ 2.55 ± 0.05, precisely 
where (within errors) /3 crosses the value 2, required for the existence 
of the first moment (normalization of the wave function). 

phase and their values are almost independent of N for the 
largest sizes explored. [ ] This is an uncommon distribution 
for the quantity x: in the Anderson model usually a = 1/2 
and the large x behavior is exponential [ ] reminiscent of 
the Porter- Thomas distribution of RMT [25]. Comparing the 
power-law tail with the exponential one of the delocalized 
phase in the Anderson problem, we conclude that already deep 
in the delocalized region, there are sign of pre-localization. 
The almost perfect collapse of the curves in the upper panel 
of Fig. 3 allows a much better finite size scaling analysis than 
any of its moments. 

As h approaches h c ~ 2.6 the elbow smoothens and a — > 1 
so that we can identify h c as the point at which a = 1, the dis- 
tribution stops being summable and necessarily the indepen- 
dence on N ceases. [26] This occurs at h c = 2.55 ± 0.05 as it 
can be seen in Fig. 4. An explicit A-dependence of <p means 
that the scaling of all the IPR's and of the diagonal entropy 
with A change abruptly and ergodicity is broken. 

The exponent (3 governs the scaling of the various IPR 9 's. 
For < q < f3 — 1 the integral in (8) is finite and IPR g ~ 
Af 1 ~ q . If instead /3 — 1 < q, since the integral in (8) is di- 
vergent the average of IPR 9 over the initial states \a) does 
not exist, but the typical value for a state should be found 
by looking at the sum of Af independent and identically dis- 
tributed variables x q a . One then finds the probability den- 
sity for J2 a <Af x a = Y (by computing and then inverting its 
Laplace transform, provided j3 > 2) as 

P(Y) oc y-^exp f-C (^Y~j j. (10) 

where j = 1 + (/3 — l)/q, (1 < 7 < 2) and C is a constant 
of O (1). This distribution has a power law tail but the typical 
value of the sum is set by the exponential as Y ^ Af 1 /^- 1 ) ^> 
Af. This implies typical values of the IPR g of a state, when 



q>P-l: 

IPRW ~ Af- q+ f^. (11) 

The different IPR g define different "critical points" h q solu- 
tions of /3(hg) = q + 1. The real transition, signaled by an 
explicit A-dependence of full distribution <fi can then be iden- 
tified by the diagonal entropy (6), or the limit as q — >• 1 of 
IPR 9 , therefore when /3 = 2. What is the possible origin 
of the power-law tail at large xl[ ] This can be linked with 
the existence of a many-body mobility edge at some energy 
E*(h), where eigenstates occupy O (Af) sites above E* and 
O (Af a ) (a(h) < 1) below E* and to a competition between 
the canonical entropy (the logarithm of the number of states 
between energy E and E + dE) and the diagonal entropy mul- 
tiplied by q. This phenomenon deserves better investigation in 
a future work. 

Summarizing, the coincident divergence of (1) (a non- 
summability of <p(x) at small x), and of (x) (non summability 
of xcj)(x) at large x) signal the beginning of the localized re- 
gion. This implies an accumulation of wave-function ampli- 
tudes towards small values typical of localized states[22]. We 
expect then that the scaling of the typical coefficients changes 
abruptly at the onset of the region in which ergodicity is bro- 
ken but the wave functions are still extended. 

This suggests a description of the localized phase in which 
a typical eigenstate is described by negligible weight on am- 
ple regions of the Hilbert space, which is reminiscent of the 
"small branching number" Bethe lattice picture of [2, 28] and 
of the eigenstates of a disordered but integrable model [29]. 



SIMILARITIES WITH AL ON THE BETHE LATTICE 

[22,30] 

Our case shows three differences from this classic topic: 
1) our lattice has connectivity O (JV) > 0(1) (but still 
<C O (Af), the volume of the system), 2) the on-site disorder 
potentials of neighboring configurations a and b are strongly 
correlated (E a — Ei = hi+i — hi <C E a , E^) and 3) our lat- 
tice is not random at all. In order to identify which of these 
three ingredients are necessary to preserve this phenomenol- 
ogy of the distribution functions we have investigated numer- 
ically a random graph with Af nodes and fixed connectivity 
N/2 and independent random energies e, on each node. We 
observe the same qualitative features in the distribution 4>(x), 
even for small h. On the contrary, for the Anderson model on 
a Bethe lattice with connectivity O (1) in the ergodic region 
we observe an exponential (or possibly stretched-exponential) 
tail at large x (the data will be presented in a future publica- 
tion). Therefore we conjecture that the necessary requirement 
for the large x power-law tail is the growing connectivity, and 
that one can get rid of the correlation of the energies and the 
specific topology of the hypercube. 

This confirms that we have the right to look at MBL as a lo- 
calization phenomenon on a Bethe lattice with asymptotically 
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large connectivity, a problem amenable of analytic treatment, 
beyond the locator expansion [28, 30]. 



SUMMARY AND CONCLUSIONS 

We have investigated the behavior of the return (or survival) 
probability as a possible detector of ergodicity breaking and of 
the MBL transition. We have shown how this question leads 
to the necessity of a thorough study of the distribution of the 
wave-function amplitudes of the eigenstates averaged over all 
energy spectrum. [ ] We then identified the major changes 
which occur to this distribution at the MBL transition point. 
The delocalized, ergodic phase looks more localized than the 
corresponding single-particle AL and RMT does not seem to 
be a good approximation, not even deep in the delocalized re- 
gion. The localized region seems very akin to the case of sin- 
gle particle AL on the Bethe lattice with connectivity O (1), 
in particular the distribution functions of the amplitudes show 
a small- a; accumulation which points towards wave functions 
localized in configuration space. We have also identified sim- 
ilarities and differences with this better studied case and sug- 
gested what are the necessary ingredients for a viable analyti- 
cal study of MBL. 
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